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y  An  experimental  and  analytical  study  was  conducted  for 
quantitative  evaluation  of  the  interdependence  of  stress, 
strain,  moisture  diffusion  and  heat  conduction  in  a  polymeric 
adhesive-bonded  layer  as  a  function  of  time.  Following 
experimental  evaluation  of  the  necessary  coefficients,  a  set 
of  nonlinear  coupled  differential  equations  is  solved 
numerically  by  the  finite  strip  method. 

It  was  shown  that  moisture  diffusion  coefficient 
correlates  well  with  the  volumetric  changes  in  the  stressed 
materials.  Experimental  data  of  the  deformation  behavior 
vs.  moisture  diffusion  in  an  assymmetrical  epoxy-aluminum 
model  exposed  to  different  hygrothermomechanical  histories 
have  shown  a  good  correlation  with  analytical  prediction.  (  ' 
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PREFACE 


The  research  on:  "Durability  of  adhesively-bonded  systems 
of  has  been  extended  over  three  years  and  comprises  of  three 
main  parts  as  follows: 

Part  A:  Hygrothermal  effect  on  the  deformational  behavior  of 
the  adhesive  material  as  a  bulk  and  as  an  in-situ 
bonded  layers . 

Part  B:  The  hygrothermal  effect  on  the  viscoelastic  and  time- 
dependent  response  of  the  bonded  system. 

Part  C:  The  third  part  which  is  presented  in  the  present 

report,  covers  the  stress-strain-time  relationship 
of  the  bonded  adhesive  as  affected  by  the  coupling 
between  the  hygrothermal  processes  and  the  mechanical 
stress-strain  field. 

The  final  goal  of  this  study  is  the  prediction  of  time- 
dependant  deformational  behavior  and  ultimate-strength  limit, 
of  the  bonded  system  under  hygrothermal  service  conditions  . 

Such  prediction  can  be  based  on  the  available  analytical 
formulation  and  experimental  findings  which  were  accomplished 
during  the  three  years  research  effort. 

The  personel  involved  in  the  research  consisted  of  three 
academic  faculty  members  working  part  time  and  two  Ph.D. 
candidates,  as  well  as  a  full  -  time  qualif  ied  high-level  technician  . 
Technical  assistance  was  provided  by  the  Laboratory  of  Mechanics 
of  Material  at  the  Department  of  Mechanical  Engineering  of 


the  Technion. 
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1 .  Introduction  1 

-  ( 

i 

The  effect  of  moisture  and  temperature  on  the  mechanical 
properties  of  bulk  polymeric  material  has  been  widely  studied. 

It  can  be  stated  generally  that  following  a  raise  in  tempera¬ 
ture  and  moisture  content,  the  mechanical  properties  usually 
undergo  accelerated  degradation. 

When  a  material  system  is  subjected  to  extra  humidity 
and  thermal  changes,  an  internal  redistribution  of  the 
moisture  content  or  temperature  develops  across  it.  Since 
the  distribution  is  nonuniform,  hvgroelastic  (thermoelastic) 
stresses  develop  and  may  persist  even  after  steady  state  is 
attained . 

It  is  a  well-known  phenomenon  that  most  polymeric 
materials  swell  due  to  moisture  absorption.  However,  when 
this  hygroelastic  expansion  is  mechanically  restrained,  the 
absorption  process  is  likely  to  proceed  differently  compared 
with  unrestrained  material. 

The  effect  of  stress  on  moisture  diffusion  can  be 
observed,  for  example,  when  imposing  a  tensile  stress  on  an 
angleplv  composite  laminate  in  a  humid  environment.  It  was 
found  by  Broutman  ei  al  [1 ,  2]  that  as  the  angle  between  the 
direction  of  load  application  ana  fiber  orientation  is 
increased  (thereby  increasing  the  tensile  stress  on  the  matrix) 
the  rate  of  moisture  diffusion  into  the  matrix  is  enhanced. 

One  of  the  first  attempts  to  evaluate  qualitatively  rhe 
effect  of  stress  on  the  moisture  diffusion  rate  was  made  by 
Fahmy  and  Hurt  [3] ;  they  developed  an  expression  for  the 
coefficient  of  moisture  diffusion  as  a  f u:  'tier,  of  the  average 
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hydrostatic  pressure  (AHP) ,  which  was  found  to  be  linear  for 
low  values  of  AHP. 

The  constraints  affecting  hygroelastic  swelling  and  the 
diffusion  process  fall  into  the  following  four  main 
categories : 

(a)  Free  swelling  -  when  the  specimen  is  very  thin  there  is 
practically  no  dependence  of  moisture  concentration  on  thick¬ 
ness  . 

(b)  Homogeneous  constraint  -  when  the  specimen  is 

thick,  a  distribution  of  moisture  will  develop  depending  on 
the  distance  from  the  outer  surface;  each  element  of  the 
material  will  tend  to  expand  according  to  its  own  moisture 
content,  with  the  attendant  effect  on  the  adjacent  elements. 

(c)  Heterogeneous  constraint  -  in  specimens  composed  of  two 
or  more  different  materials,  (i.e.  bonded  systems,  composites 
etc.)  the  restrained  hygroelastic  swelling  has  two  main  causes 

1 .  Moisture  differentials  due  to  different  dif f usivities 
and  2.  Difference  in  hygroelastic  coefficients  of  expansion. 

(d)  External  application  of  mechanical  and/or  thermal  loading 

The  experimental  part  of  the  present  research  deals  with 
the  effects  of  the  above  constraints  on  the  diffusion  process. 
The  analytical  study  is  mainly  concerned  with  the  interactions 
between  the  three  different  processes  and  mechanisms  which  are 
active  simultaneously,  namely;  heat  conduction,  mechanical 
response  and  moisture  diffusion.  Several  studies  evaluated 
the  extent  to  which  these  processes  affect  each  other. 

Sih  and  Shih  [3,4,5]  take  account  of  the  combined  effects 


of  moisture  and  temperature  by  expressing  the  moisture  diffu¬ 
sion  coefficient,  D,  as  an  exponential  function  of  the  ab¬ 
solute  temperature,  and  by  inserting  coefficients  linking 
moisture  and  temperature  into  the  diffusion  equations.  The 
effect  of  temperature  on  the  hygroelastic  expansion  coeffi¬ 
cient,  8,  (HEC)  and  the  effect  of  moisture  on  the  thermo¬ 
elastic  expansion  coefficient,  a,  (TEC)  has  been  evaluated 
experimentally  by  Ishai,  Gali  and  Yaniv  [6]. 

The  combined  influence  of  stress  and  temperature  is 
expressed  through  considering  the  rate  of  plastic  work  as  a 
distributed  heat  generation  source.  In  case  of  slow  rates  of 
heat  transfer  from  the  stressed  body  to  the  environment,  the 
temperature  is  likely  to  rise.  This  in  turn  may  affect  the 
mechanical  properties  and  therefore  the  state  of  stress  [7,8]. 

Chung  and  Bradshaw  [9]  have  developed  a  set  of  inter¬ 
dependent  equations  governing  moisture  diffusion,  heat  con¬ 
duction,  and  mechanical  behavior.  In  order  to  describe  the 
mechanical  behavior  of  polymers,  a  kernal  function  of  a 
viscoelastic  nature  was  adopted.  Further  experimental  data 
and  analytical  studies  on  this  subject  have  been  reported  by 
Sprinqer  et  al  [10,11]. 

When  a  polymeric  material  is  stressed  as  a  phase  in  a 
bonded  system  or  a  composite  material,  it  may  exhibit  a 
pronounced  time-dependent  elasto-viscoplastic  behavior.  Such 
behavior  is  usually  sensitive  to  hyqrothermal  effects.  Due 
to  the  complexity  of  the  processes  involved,  there  seems  to 
be  no  solution  available  to  date  which  rigorously  considers 
the  interdependence  of  ana  the  links  between,  the  time- 
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dependent  mechanical  behavior  and  the  hygrothermal  processes. 
The  present  work  is  an  attempt  to  fill  this  gap  by  tackling 
the  subject  both  theoretically  and  experimentally. 


2.  Experimental  Investigation 

2.1-  Effect  of  external  load  on  moisture  diffusion 

The  effect  of  stress  on  moisture  diffusion  can  be  deter¬ 
mined  by  measuring  and  calculating  the  amount  of  water  gain 
per  unit  weight  of  the  material  in  its  dry  state,  as  a 
function  of  the  loading  level  (see  Appendix  A) . 

2.1.1:  compressive  loading 

Two  kinds  of  polymers  were  examined  -  Epoxy  resin 

® 

(Epon  828  and  Varsamid) ,  and  FM-73  ,  which  usually  serves  as 
structural  adhesive  film.  Two  types  of  specimens  were 
machined  from  cured  bulks  -  thin-walled  tubes  for  low  and 
moderate  loads  (up  to  8.0  MPa)  and  thick-walled  ones  for  high 
load  (up  to  25.0  MPa),  as  shown  in  Fig.1. 

The  level  of  the  externally  applied  compressive  load  was 
determined  by  presetting  the  length  of  the  spring  (1)  in  the 


device  shown  in  Fig. 2.  By  releasing  the  nut  (4) ,  the  force 
exerted  by  the  spring  is  transferred  through  the  arm  (5)  and 
disk  (2)  to  the  specimens  (3) .  The  spring  coefficent  was 
designed  to  maintain  the  applied  load  within  less  than  5% 
change  due  to  swelling  and  creeping  of  the  specimens  during 
the  test.  The  device  can  accommodate  three  to  four  specimens 
in  series,  a  fact  that  ensures  an  identical  force  on  all  of 
them.  In  order  to  minimize  interruptions  of  the  test  due  to 
readings,  rapid  removal  of  the  specimens  is  effected  by  closing 
nut  (4)  and  opening  nut  (6) . 

Six  devices  under  different  loads,  comprising  specimens 
preliminarily  dried  in  a  vacuum  oven  at  80°C  were  assembled 
in  a  distilled  water  container,  which  was  kept  in  an  oven  at 
40°  ±5°C  (Fig. 3)  for  different  durations.  Weight  and  dimen¬ 
sional  change  readings  were  taken  at  regular  intervals,  and 
compared  to  a  reference  set  of  specimens  kept  dry  at  40°C. 

An  additional  set  was  kept  undisturbed  in  the  water  container 
throughout  the  test  in  order  to  examine  the  effect  of  removal 
of  the  specimens  on  moisture  uptake.  Test  conditions  for  all 
variants  are  summarized  in  Tables  1  through  3. 

Figs.  4  to  6  show  the  relative  moisture  absorption  vs. 
exposure  time  at  the  different  load  levels.  It  is  clearly 
seen  that  the  higher  the  compressive  load  level,  the  lower 
the  moisture  absorption.  Because  of  the  relative  shortness 
of  the  test  the  long-term  effect  of  compressive  load  on  the 
maximum  amount  of  water  gain  cannot  be  conclusively  estab¬ 
lished,  but  a  general  picture  can  be  derived,  as  outlined  in 


Section  3. 


2.1.2:  tensile  loading 


The  tensile  specimens  were  machined  to  a  dog-bone  shape 
from  a  cured  FM-73  plate.  Four  different  lenaths  of  the  inner 
zone  (I)  were  used  (Fig. 7)  in  order  to  isolate  the  deformation 
in  it  from  the  nonuniform  stress  field  at  the  end  zones  (E) , 
which  were  identical  in  all  specimens.  Each  device  accom¬ 
modated  four  specimens.  The  tensile  test  procedure  was 
similar  to  its  compressive  counterpart  using  the  device  shown 
in  Fig. 8.  Test  conditions  are  summarized  in  Table  4.  Moisture 
content  vs.  exposure  time  at  different  load  levels  is  shown 
in  Fig. 9.  The  procedure  for  determination  of  moisture  content 
in  specimens  under  tension  is  described  in  Appendix  A. 

2.2-  The  effect  of  homogeneous  material  constraint  on 
diffusion  of  moisture 

For  this  series  five  sets  of  specimens  were  prepared 
from  cured  FM-73  plate  of  different  thickness.  The  specimens, 
also  preliminarily  dried  in  a  vacuum  oven  at  80 °C  for  seven 
days,  were  tested  with  periodical  examination  as  under  (2.1). 
Test  duration  was  14  days. 

Fig. 10  shows  the  relative  moisture  content  vs.  exposure 
time.  A  significantly  higher  moisture  absorption  rate  was 
found  for  the  thinnest  specimen,  which  will  serve  later  as 
constraint-free  reference. 

2.3-  The  effect  of  heterogeneous  restraint  on  moisture 
diffusion 

The  second  objective  of  the  experimental  phase  was  to 
obtain  experimental  data  in  which  the  effect  of  the  reciprocal 
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influences  between  moisture,  temperature,  and  mechanical 
behavior,  will  be  represented.  These  data  were  used  to 
verify  the  analytical  results. 

Three  types  of  asymmetrical  doubler  (NSD)  specimens  were 
manufactured  (Fig. 11),  the  material  composition  and  dimensions 
of  which  are  listed  in  Table  5.  In  specimens  of  series  C-1 
and  C-2,  the  FM-73  polymeric  layer  was  cured  at  123°C  during 
specimen  manufacture,  while  in  series  T- 1 ,  the  polymeric 
layer  was  first  cured  separately  and  was  later  cemented  to 
the  aluminum  layer  and  cured  at  75°C.  After  curing,  the 
specimens  had  become  warped  (Fig. 11)  due  to  the  difference 
in  the  TECs  of  the  material  constituents .  Because  of  the 
higher  thermal  loading  the  central  deflection,  and  thus  the 
compressive  stresses  in  the  outer  strata  of  the  polymeric 
layer  of  C-1  and  C-2  specimens  were  greater  than  those  in  the 
T-1  series. 

All  specimens  were  immersed  in  a  water  container,  which 
was  placed  in  an  oven.  At  certain  time  intervals  the  central 
deflection  (Fig. 11)  and  water  uptake  was  measured.  Due  to  the 
diffusion  of  moisture  the  polymeric  layer  swelled,  a  fact  that 
was  reflected  as  a  decrease  in  the  central  deflection.  The 
experimental  results  shown  in  Figs.  12  to  15,  clearly  show 
that  the  rate  of  diffusion  in  the  specimens  with  the  less 
compressed  outer  layer,  i.e.  series  T-1,  is  higher  than  in 
series  C-1  and  C-2.  Furthermore,  between  the  11th  and  the 
13th  day  of  exposure,  a  deviation  from  the  initial  slope  of 
the  T-1  curve  in  Fig. 13  is  evident,  caused  by  an  increase  in 
the  value  of  the  moisture  diffusion  coefficient.  Examination 


f 


of  the  deflection  vs.  time  curve  (Fig. 14)  reveals  that  at 
the  stage  referred  to  the  curvature  of  T-1  specimens  changes 
direction.  As  a  result  of  the  bending,  the  stresses  on  the 
outer  side  become  tensile  instead  of  compressive,  but  no  such 
reversal  occurs  in  series  C-1  and  C-2.  This  explains  the 
increased  rate  of  moisture  diffusion  in  T-1  specimens  shown 
in  Figs.  13  to  15. 


3 .  Analytical  Formulation  For  the  Moisture  Process 


Isothermal  moisture  diffusion  in  polymers  can  be 
described  by  Pick's  lav;  [12J,  as  follows: 


D 


x 


8C 

at 


(3.1) 


with  boundary  and  initial  conditions: 


where  D 

x 


C 

c 

c, 

h 


0  <  x  <  h  t  £  0 

x  =  0  ,  h  t  >  0 

diffusion  coefficient  in-x  direction, 
moisture  concentration 
initial  concentration 
final  concentration  (saturation) 
plate  thickness 


(3.2) 


Eq.(3.1),  solved  subject  to  condition  (3.2),  yields  the 


well-known  expression: 


C-C  . 
1 


(2  j  +  1  )  TTX 


exp  - 


(2 j  +  1 )  Tt  D  t 


(3.3) 


For  a  given  moisture  distribution,  the  water  uptake  (in 
weicrht)  is  qiven  by: 


ir.  =  nS.  ,  C  dx 


(3.4) 


Taking  into  account  the  finite  dimensions  of  the  plate 
(n  -  width,  l  -  length),  integration  eq.(3.3)  yields: 


m  -m . 
m  i 


—  =  1  -  — ?  1  (2j  +  1)  exp  —  ( 2 j  + 1 )  ^tt 

_ “  • 


TT  j  =  0 


2_2  Dt 


,  (3.5) 


where 


D  -  Dx(1  *  *  ♦  f)2 


(3.6) 


The  total  weight  of  the  specimen  will  be  then: 


W  =  +  m 


(3.7) 


and  the  relative  water  gain,  M,  is  given  by: 


“d  _  _m_ 
Wd  Wd 


(3.8) 


where  is  the  dry- state  weight. 

For  short  exposure  periods*,  where  ^  <  1,  using  eq.(3.8) 

h 

equation  (3.5)  can  be  approximated  by: 


*  A  typical  short  period,  for  h  s  2mm  plate,  is  about  10  days 
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Figures  16  to  18  show  the  RMC  vs  /t  at  different 
compressive  loads,  bearing  out  the  predictions  of  eq.(3.9). 
The  higher  the  compressive  load,  the  lower  the  diffusion 
rate,  i.e.  the  CMD. 

The  plots  in  figs. 16-1 8  permit  determination  of  the 
ratio  >l(a)/M(0)  and  the  nondimensional  diffusivity  D(cj)/Dq 
can  be  calculated  eq.(3.11)  and  plotted  vs  stress.  Fig. 19 
indicates  that  at  low  compressive  stresses*  there  is  a  linear 
relation  between  CMD  and  stress  which  can  be  expressed  as: 

D(o)  =  Dq ( 1 +kJ^ )  .  (3.12) 

being  the  first  invariant  of  the  general  stress  tensor  and 
k  a  constant.  The  present  values  for  eq.(3.12)  are: 

D (a)  =  D0 [ 1+ (9. 1 2  ±  ,32)1 Q~2  J  .  (3.13) 

The  effect  of  tensile  load  on  the  increase  in  diffusion 
rates  is  shown  in  fig. 20.  Plotting  the  relative  diffusivity 
vs  stress  show  that  the  linear  relation  again  exists  only  at 
low  stress  levels  (fig. 21).  The  nonlinearity  at  higher 
stresses  is  attributable  to  the  visco  plastic  mechanical 
behavior  rather  than  to  stress-diffusion  coupling  as  will  be 
discussed  later. 

If  the  "Relative  Reduced  Volume"  (RRV)  is  defined  as 
the  difference  between  the  relative  total  volume  change  and 
the  free  hygroelastic  swelling  per  unit  volume,  the  follow¬ 
ing  may  be  stated:  The  greater  the  RRV,  the  greater  the  rate 
of  moisture  diffusion,  and  vice  versa;  in  other  words,  if  the 

*  But  not  at  high  ones  -  see  [3]  and  [15] . 
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=  2 
M  /i  h 
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Based  on  the  above  assumptions,  a  linear  relation  between  the 

relative  moisture  content,  M,  and  the  square  root  of  the 

coefficient  of  moisture  diffusion  (CMD) ,  is  obtained.  The 
/t 

plot  of  M  vs  -r-  for  short  periods  should  be  a  straight  line 
the  slope  of  which  is  directly  proportional  to  /d  . 

The  ratio  of  the  relative  moisture  content  (RMC)  of  a 
stressed  and  an  unstressed  specimens  (of  equal  thickness)  at 
a  given  time  is  as  follows: 


M(o)  _  Mm(o)  |  D  ( g ) 
M(0)  M  (o)  Dn 


(3.10) 


where  (o)  and  (0)  refer  to  the  stressed  and  unstressed 
conditions,  respectively. 

In  earlier  studies,  [13]  and  [14],  it  was  found  that 
external  loading  usually  does  not  affect  the  saturation 
content,  but  the  diffusion  rates  only.  Since  loads,  in  the 
present  experiments  are  relatively  low,  this  assumption  can 
ue  employed  here  to  give: 


M (a)  D(g) 
M(0)  '  D. 


(3.11) 
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total  state  of  stress  causes  a  decrease  In  the  RRV,  the  rate 
of  diffusion  may  be  expected  to  decrease.  The  experimental 
results  (Fia.21)  and  the  fact  that  in  a  state  of  yielding  or 
plastic  flow  the  change  in  the  relative  volume  usually 
approaches  zero,  confirm  the  above  assumption  about  the 
dependence  of  the  diffusivity  on  RRV.  This  justifies  the 
expression  of  the  diffusivity  in  terms  of  strains  rather 
than  by  relating  it  directly  to  stresses.  Such  an  approach 
permits  a  simple  formulation  of  the  link  between  stresses 
and  moisture  diffusion,  as  will  be  shown  later  on  . 

The  effect  of  homogeneous  constraint  on  moisture 
diffusion  is  shown  in  ficr.22.  In  spite  of  the  large  differ¬ 
ences  in  specimen  thickness,  diffusion  rates  are  almost  the 
same.  It  can  be  assumed,  then,  that  the  developed  low 
hygroelastic  stresses  have  a  negligible  effect  on  moisture 
diffusivity.  An  equivalent  hygroelastic  stress  was  computed 
(see  Appendix  B) ,  and  the  results  are  shown  in  Fig. 23.  The 
two  figures  indicate  that,  for  most  practical  purposes,  the 
effect  of  homogeneous  material  constraint  on  the  CMD  is 
negligible . 

Better  insight  into  the  nature  of  the  mutual  influence 
of  stress  and  of  moisture  diffusion  may  be  achieved  by 
examinina  some  of  the  results  in  the  micro  level.  After  the 
empty  spaces  of  the  specimen  have  been  filled  with  water, 
further  intake  is  conditional  on  its  expansion,  which  in  turn 
entails  overcoming  the  intermolecular  forces.  Resistance  of 
the  material  to  swelling  which  increases  with  the  amount  of 


pm*: 
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additional  moisture  diffusion  is  effected  by  application  of 
a  compressive  load.  By  contrast,  tensile  load  tends  to 
release  the  resistance  to  swelling  and  the  diffusion  rate 
increases . 

The  difference  between  the  total  relative  volume  change 
and  the  free  hygroelastic  swelling  per  unit  volume  was 
defined  as  the  Relative  Reduced  Volume  (RRV) .  As  the  net 
result  of  mechanical  and  hygrothermal  loading  is  reduction 
or  increase  of  the  RRV  of  a  certain  unit  volume,  so  the  rate 
of  moisture  diffusion  is  also  reduced  or  increased  respec¬ 
tively.  This  postulation  is  in  agreement  with  the  nonlinear 
stress-diffusion  relation  mentioned  earlier.  It  is  well 
known  that  the  volumetric  changes  in  a  material  under  yield¬ 
ing  or  flow  usually  tend  to  zero  and  so  do  the  changes  in  the 
RRV.  Thus  the  diffusivity  approaches  a  constant  value  at 
high  stress  levels. 

4.  Mathematical  Formulation  For  the  Hyarothermal-Elasto- 


Viscoplastic  CouDling  Effect 


The  differential  equations  from  which  moisture  and 
temperature  distributions  are  determined  can  be  obtained 
from  the  basic  laws  of  mass  conservation  and  energy  conser¬ 
vation  (first  law  of  thermodynamics)  and  from  the  inequality 
of  the  second  law  of  thermodynamics  [16].  The  equations 
necessary  for  describing  the  mechanical  response  can  be 
obtained  from  the  law  of  momentum  conservation,  and  by 
assuming  a  certain  model  of  viscoplastic  behavior.  These 


basic  laws  lead  to  the  well  known  heat  conduction,  moisture 
diffusion  and  equilibrium  equations. 


4.1  -  Diffusion  equations 


The  linked  form  of  the  heat  conduction  and  moisture 
diffusion  equations  in  a  2-D  rectangular  domain  (Fig. 24) 
is  expressed  as: 


S'c'Vti’  *  ^(CDyff>  -  ^IC'XTT>  *  0 

■— (  D  )  +  -dL(  D  — ^ )  -  JL-(t-X  C-  —  VI 1  -  0 
3x  T  x3x'  3y  T  y3y  3t(T  Ac"  pcW)  0 


(4.1) 
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where  TD^  is  the  thermal  diffusivity  and  CD^  is  the  coefficient 
of  moisture  diffusion  in  the  i-th  direction  (i  =  x,y):  x„  anc 
XT  are  coefficients  linking  moisture  and  temperature;  o  and  c 
are  respectively  the  density  and  the  specific  heat,  and  W  is 
the  heat  generated  within  the  domain. 

The  initial  conditions  are: 

U (x ,y , 0)  =  0 (x,y)  (4.2) 


where  U  represents  either  T  or  C. 

The  boundary  conditions  can  be,  for  example: 
U  (x ,  a ,  t)  =  UQ2  on  r2 


and  due  to  symmetry 


du  = 

3y 


0 


on  f  ^ 


(4.3) 


(4.4) 


where  UQ2  is  a  given  value  on  r2  . 

The  dimens ionaless  variables  are  defined  as 

D 

,  .  X  _  y  _  R  , 

1  =  k  •  n  =  t  ;  T  =  ri  1 

n 


where  t  is  dimensionaless  time  and  D  is  a 

R 

value  of  diffusivity.  The  physical  variabl 
shifted  and  scaled  as  follows: 


(4.5) 

certain  reference 
es  C  and  T  are 
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U  -  U 


U  = 


02 


Ur  -  U 


02 


where 


UD  =  max [U (x) ] 
K  xefl 


(4.6) 


(4.7) 


No  special  symbols  are  used  to  distinguish  dimensional 
variables  from  dimensionless  ones.  The  context  affords 
sufficient  means  to  avoid  confusion.  Thus,  the  dimensionless 
form  of  the  diffusion  equations  is  given  by: 


tc<cdx  +  ^2^(cdyS)  -  ^(c  ■  xtt)  =  0 


3C, 


9  t  n  11)  +  il)2— (  D  — )  -  — (T  -  X  C-XwW)  =  0  (4.8) 
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The  uniform  initial  conditions  are: 
UK,  n,  0)  =  1 


(4.9) 

(4.10) 


and  boundary  conditions: 

U(0,  n,  T)  =  u(i ,  n,  T)  =  UK,  i ,  t)  =  o 

on  r2,  r3,  r4  (4. id 


and 


9U 

9n 


o 


on  r1 


(4.12) 


When  moisture  and  temperature  fields  are  approximately 


described  by  C  and  T  respectively,  the  coupled  nonlinear 
diffusion  equations  can  be  transformed  into  the  weak  form 
by  Galerkin's  formulation,  as  follows: 


3(G  )  h  ar  2 

.0..  I?  *  *2 


Q 


3£  c  x  3£ 


3  (G  ) 
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3n  c  y  3n  3t  t  c  i 
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where  (Gc)i  and  (GT) ^  are  the  i-th  test  function  of  C  and  T 
respectively. 


4.2-  The  plasto- viscoplastic  model 
4.2.1:  basic  equations 


The  general  strain  tensor,  c,  can  be  expressed  as: 

( 4. 14) 


=  +  !vP  +  5t  +  £h 


where  the  constituents  are  the  elastic,  viscoplastic,  free 
therinoelastic ,  and  free  hygroelastic ,  strain  tensors  respects 
vely.  The  elastic  strain  tensor  is  related  to  the  general 
stress  tensor  o,  by: 


2  -=  ?£e 


(4.15) 


where  D  is  the  matrix  of  elasticity  which,  in  the  general 
case,  is  moisture  and  temperature  dependent.  Assuming  that 
8  and  a  are  independent  of  moisture  and  temperature  respecti¬ 
vely,  in  the  specific  ranges,  the  hygrothermal  strain  tensors 
can  be  written  as: 

5h  '  ?6C 


e_  =  aAT 
-T 


(4.16) 
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Viscoplastic  flow  occurs  when  a  scalar  yield  function, 
F,  is  positive  (Owen  and  Hinton  [18]) : 

P  -  F(0,evp.  T,  C)  -  P0  >  0  (4.17) 

where  F  is  the  equivalent  stress  in  accordance  with  the 
accepted  yield  criterion,  and  F0  is  ike  jjrield  stress  which, 
in  the  case  of  linear  strain-hardening,  Has  the  form: 

Fo  =  °Y  +  ^vp  .  (4.18) 

Here  k  is  the  strain-hardening  parameter,  is  the  initial 
yield  stress,  and  is  the  effective  viscoplastic  strain  - 

defined  by: 

Svp  ;  ‘Kp’ij'Vij1*  <i,5  =  1,2,3)  (4.19) 


Here  and  in  the  following  the  summation  convention  is  used. 

As  in  plasticity,  the  flow  law  of  viscoplasticity  can 
be  defined  as: 


The  dot  (’)  indicates  time  differentiation;  Q  is  the  plastic 
potential,  which  has  the  value  of  the  equivalent  stress,  F. 
This  approach  is  similar  to  that  used  in  the  associated 
theory  of  plarticity.  The  viscoelastic  multiplyer,  X,  is 
given  by: 

X  =  y  * (F)  (4.21) 


where 

C  F  >  0 

♦(F)  =.  (4.22) 

0  F  S  0 

m 


and  y  is  the  fluidity  parameter.  The  following  viscoplastic 
flow  function  is  used: 


*  »  (F/Fd)N  .  (4.23) 

The  model  parameters,  y  and  N,  are  determined  experimentally. 
With  the  above  assumptions,  the  viscoplastic  strain-rate 
tensor  becomes: 

e  =y  MF)  (4.24) 

~VP  oC 


It  depends  on  the  equivalent  stress  F,  whose  value  is  deter¬ 
mined  depending  on  what  yield  criteria  are  employed.  Most  of 

the  commonly  used  criteria  can  be  expressed  thus: 

F  *  c1J1  +  c2{J2)i  (4.25) 

where  J1  is  the  first  invariant  of  the  general  stress  tensor, 

I 

J  2  is  the  second  invariant  of  the  deviatoric  stress  tensor, 
and  c^  and  c2  are  constants  which  define  the  specific  yield 
criterion. 


Rather  than  the  classical  Von-Mises  and  Drucker-Prager 
criteria,  a  modified  yield  criterion  is  used  in  this  work  t6i 
This  criterion  is  suitable  for  specially  ductile  materials 
such  as  the  FM-73  employed  here.  The  constants  in  this  case 


are : 


/3(X  +1) 
m 


(4  .26) 


\  being  the  ratio  between  the  tensile  and  the  compressive 
xn 

strengths . 
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The  principle  of  virtual  work  in  the  absence  of  body 


forces  is  written  as: 


/Se^odfi  -  /  6u^fdr  =  0 


(4.27) 


where  u  =  {u,v}  is  the  displacement  vector,  f  is  an  external 
traction  vector,  and  is  the  boundary  portion  of  the  domain 
Q,  the  tractions  on  which  are  given.  Using  the  strain  dis¬ 
placement  relation  e  =  Bu  produces  the  equilibrium  equations: 


/  B  crdft  -  F  =  0 
~ 

in  which  the  external  force  vector  is  given  by: 

F  a  /  fdT 
~  Tr~ 


(  4.28) 


( 4.29) 


In  the  general  case  the  matrix  operator  B  is  composed  of  a 
linear  and  a  nonlinear  part,  but  since  the  role  of  the  non¬ 
linear  part  is  negligible  when  small  displacements  are 
considered,  it  has  been  omitted  in  this  work. 

The  incremental  form  of  equilibrium  equations,  for  the 

n-th  time  interval,  At  ,  is  written  as: 

n 


/BTAondn  -  AFn  =  0 


(4.30) 


where  Ao11  and  AFn  are  the  changes  in,  respectively,  the  stres 
tensor  and  in  the  external  load  vector,  during  the  time  inter 
val  A t ^ ,  which  is  defined  by: 


Y 


i 
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The  relation  between  stress  and  elastic  strain  increment  is: 


Aan  =  DnAcn 
-  -  ~e 


or,  using  equation  4.14: 


Aon  =  Dn[Aen  -  Ae^p  -  Ae£  -  Aej|]  . 


(4.32) 


(4.33) 


Substitution  of  equation  4.33  into  equation  4.30  and  applying 
the  strain  displacement  relation,  yields  the  following  equili¬ 
brium  equations: 

/BTDnBAundfi  -  (AFn  +  /BTDn(Ae”  +  Aejj-Jdfl]  =  0  (4.34) 

Si"  ~  ~  ~VP  ~HT 

where 


**St  *  4sS  *  *?;  • 


(4.35) 


4.2.3;  the  mutual  influence  of  moisture  diffusion  heat 
conduction,  and  mechanical  behavior 

Moisture  and  temperature  affect  both  the  state  of  stress 
and  deformation,  directly  through  hygrothermal  strains  and 
indirectly  by  affecting  mechanical  properties.  However, 
the  mechanical  behavior  may  also  affect  the  moisture  and 
temperature  distributions  within  the  system. 


During  the  viscoplastic  process,  heat  is  generated 
through  dissipation  and  is  assumed  to  be  proportional  to 
the  rate  of  the  specific  viscoplastic  work,  W^p  l 7] . 

While  the  effect  of  this  heat  generation  source  can  be 
neglected  in  metals  due  to  their  high  thermal  conductivity, 
in  polymers,  in  which  thermal  conductivity  is  relatively  low. 
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temperature  profiles  may  be  noticeably  affected  by  this  source. 


In  section  3  it  was  stated  that  the  tensor  of  the  moisture 

(  , 

diffusion  coefficients  D  is  affected  by  the  changes  in  the 
Relative  Reduced  Volume  (RRV) .  Let,  then,  the  reduced  strain 
tensor  be  defined  as: 


e  =  e  -  e 


H 


=  e  +  e  +  cm 
~e  ^vp  _T 


( 4 . 36  ) 


hence,  the  RRV  can  be  expressed  as: 


V 


R 


e  .  . 
ll 


(4.37) 


In  the  present  work  the  discussion  has  been  restricted 
to  isotropic  materials,  and  the  only  nonzero  terms  of  cD  are 
therefore  the  diagonal  terms,  which  are  described  by: 

cDii  =  (cDo,iit1  +  xo(^R,i1  (no  sum)  (4.38) 


where  (  D  ) . .  is  the  i-th  diagonal  term  of  D  in  the  stress- 
free  reference  state,  X  is  a  constant,  and  (S  ) .  is  defined 

D  K  1 

by: 


<Vi 


-=  VR  - 


ii 


(no  sum) 


(4.39) 


The  physical  meaning  of  (SR)i  is  the  relative  change  in  the 

area,  the  normal  of  which  is  colinear  with  the  direction  of 

the  i-th  diffusion  component.  In  a  like  manner  as  RRV,  (S  )  . 

R  X 

can  be  called  the  Relative  Reduced  Surface  -  RRS . 


To  sum  up,  the  governing  equations  to  be  solved  simul¬ 
taneously  are  the  diffusion  equations: 


C  1  ' 


3  (G_) 


rr  V  i  9T  2  3  (CV  i  »T  ?  ' 

— - *TDy^I7IT->cc-1pwvp»'GT»i)dB  ■  0 


and  the  equilibrium  equations: 


/BTDnBAundfi  -  [AFn  +  /BTDn(Acfr  +  Ae[L)dfi]  =0  .  (4.41; 

Q'  *  '  *  -  A"  *  *  P  'HT 

\ 

The  above  set  of  equations  is  coupled  and  nonlinear  as  regards 
both  time  and  space.  The  numerical  solution  of  these  equations 
will  be  described  in  the  next  section. 

5 .  Computational  Model 


The  physical  problem  of  the  present  model  is  composed 

of  three  different  but  simultaneous  processes,  viz.  heat 

conduction,  mechanical  response,  and  moisture  diffusion. 

Each  process  has  a  typical  time  scale  -  the  temperature  time 

_2 

scale  is  of  0(10  min),  viscoplastic  response  is  of  0(1  min) 
and  diffusion  time  scale  is  of  0(10^  min).  Hence,  time  step 
size  is  automatically  changed  in  accordance  with  the  respective 
time  scales.  For  the  convenience  of  the  numerical  solution,  th 
diffusion  equations  were  solved  first  at  each  time  step, 
followed  by  the  equilibrium  equations.  The  flow  charts 
describing  the  computational  procedure  are  given  in  Appendix  C. 

5.1-  Diffusion  equations 


According  to  the  Finite  Strip  Method  (FSM) ,  the  repre¬ 
sentations  of 'moisture  and  temperature  distributions  within 


the  given  domain  (Fig. 25)  are  approximated  by: 

,  m 


U  = 


N±  ( e )  Ym  ( n ) 


(5.1) 


where  and  Y  are,  respectively,  the  local  and  the  global 


trial  functions;  Cm  and  T“‘  are  the  approximated  values  of  C 


„m 


and  T  at  the  i-th  nodal  point,  corresponding  to  the  m-th 
global  trial  function.  The  summation  in  equation  5.1  is 
over  all  I  nodal  points  and  M  global  trial  functions. 
Substitution  of  the  above  approximation  into  the  Bubnov- 
Galerkin  formulation  of  eqs.4.41  yields  the  following 
equations  for  the  i-th  nodal  point: 


1  3  (G  )  i  ir  t  3  (G  )  .  o  '  ~  : 

'  '  -aPcDJ  *  *  aW„  =  o 

rpO  C  C  X 


c  y 


C  1 


1  3  (G_)  .  -sm  „  3  (G_)  .  •*  «  - 

/  f  — rr-i  D  I?  +  ^ D  ~  (T-x„C-a.w„J (cm) ..  c^dn 


,  3£  Ty3n  3t 


p  vp  T  i 

(5.2) 


Essential  boundary  conditions  are: 
UU,1  ,t)  ■  0  on  r2 


(5.3) 


Since  this  condition  must  be  satisfied  for  every  £  ,  it  is 
required  that: 

m  ,  c  . . 

Y (1 )  «  0  (m  =  1  ,2, . . . ,M)  (5.4) 


The  following  set  of  complete  functions  satisfies  this 
requirement : 


Ym(n)  =cos(umn) 


(  5.  5) 


where 


2m- 1 
ym  =  “T" 


(m  =  1  ,  2  ,  .  . .  ,  M) 


(5.6) 


Boundary  conditions  of  symmetry  are  also  fully  satisfied 
exactly  by  this  approximation,  although  this  is  not  strictly 


necessary , 


Substituting  approximation  (5.1)  into  (5.2)  gives  the 


following  semidiscrete  equation: 
KU  +  CU  =  F 


(5.8) 


where 


u 


1  T  2  T 
=  nu1)  ,  (UZ)  , 


T  T 

,  (U1 )  ,  (U*)  ] 


<  T 


and 


A  * 


(U1) 

=  [C1  ,  T  ^  ,  C2  > 

T  2  /  • 

11  12 
/  A  , 

A1  3 

000$  **  /  •  0  0  f 

a1m 

21  .22 
•  A  , 

A2j 

000$  /  0  0  0  f 

a2m 

il  ni2 

i  A  , 

0  •  0  f  /  0  0  0  f 

A1M 

Ml  *M2 
k.  *  A  , 

AMj 

0  0  0  $  **  $  0  0  0  $ 

amm 
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Here  A  represents  the  structure  of  both  the  K  <_nd  the  C 
coefficient  matrices. 

According  to  Hughes'  implicit  time  marching  scheme  117], 
equation  4.8  can  be  written  at  time  Tn+e  =  (n+6)Ax  as: 


-n+e-n+e^n+eVn+e  =  ?n+e 


(5.10) 


where : 


<D 

VI 

O 

<  1 

Vn+e 

*  ZT'Sn.l 

-  V 

Ve 

*  n-9lVn 

+  eVl 

^n  +  e 

-  (,-e,?n 

+  e?n  +  1 

*n+e 

Tn  +  0  ^ 

c  . 

=  C(U_._, 

T _ ) 

(5.11) 
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By  substituting  eq.  5.11  into  eq.  5.10  and  after  further 
manipulations,  a  set  of  algebraic  equations  for  time  step 
n+1  is  produced: 


1 0K 


n+e 


J_  c  1 

At  ~n+9  J 


Vi 


=  I ( e- i ) k 

~n+e 


IT  W  Vn  *  ? 


n+  6 


(5.12) 


These  nonlinear  equations  are  solved  by  a  semi-iterative 
technique  based  on  the  Newton-Raphson  method.  At  each 
iterative  stage,  the  linearized  equations  were  solved 
directly . 

5.2-  Equilibrium  equations 

The  incremental  form  of  the  equilibrium  equations  was 
expressed  in  this  work  as: 

/BTDnBAUndn  -  [AFn  +  /  bViAe"  +  AeJLjdtl]  *  0  (4.34) 

Q  n  vp  -Hi 

in  which  the  hygrothermal  strain  increment  tensor  is  given 
by: 

Ae£t  =  §ACn  +  aATn  .  (5.13) 

The  incremental  viscoplastic  strain  tensor,  Ae11  ,  is  a 

~vp 

function  of  the  stress  itself;  hence  the  form  of  equation 
4.34  is  implicit  and  has  to  be  formulated  numerically 

for  the  solution  of  this  equation.  Using  an  implicit  time 
stepping  form  (16],  a is  expressed  as: 

Acn  =  AtJ  ( 1  -6 )  en  +  6en+1]  .  (5.14) 

'vp  n  ~vp  ~vp 

A  Taylor  series  expansion  is  usefully  employed  to  evaluate 


/B  D  Bdn 
n~ 

/BTDn(Ae”T  +  £vp^"n)d«  +  n 


(5.23) 


In  order  to  evaluate  D  ,  the  pseudo-viscoplastic 
compliance,  Hn,  which  is  defined  as: 


„  3  e“ 

tn  _  -vd 


(5.24) 


must  be  determined  explicitly. 

Let  the  flow  vector,  a,  be  defined  as  ([181): 


(5.25) 


then  H  can  be  expressed  in  terms  of  a,  using  the  chain  rule 
of  differentiation  and  the  definition  of  ,  by: 


h  -  i  *13  ^  3  $  T, 

H  -  Y  { }  . 


(5.26) 


As  already  mentioned,  yield  criteria  in  this  work  are 
assumed  to  have  the  general  form:  F  =  c 

appropriate  differentiation  of  F  with  respect  to  o,  and 
after  some  simple  mathematical  manipulations  one  gets: 


H  =  p.Mi  (i  =  1,2, 3, 4) 


in  which,  for  plane  stress  problems 


(  5.27) 
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where  S  ,  S  etc.  denote  the  deviatoric  stress  tensor 
xx  yy 

components ,  and 

C2  2 

p1  =  yl~(J'2r  *  } 

4  c 

p2  =  [c2||  - 

,2  3$  . 

P3  '  Y  {  1  3F  } 

p4  =  y{c1c2(J^)”i  •—  )  (5.29) 

A  spatial  approximation  of  the  equilibrium  equation  is 
obtained  by  finite-strip  formulation  of  the  displacements  u 
and  v,  in  the  x  and  y  directions  respectively,  as  follows: 


u  -  u  =  u^(T)Ni(OVm(n) 

v  ~  v  =  v™  (  t  )  ( £ )  Zm  ( n )  (5.30) 

where  are  local  trial  functions,  Y^Cn)  and  Zm(n)  are 
global  trial  functions,  and  u™  and  v™  are  the  aporoximated 
values  at  the  i-th  nodal  point,  corresponding  to  the  m-tn 
global  trial  function.  The  global  trial  functions  are  given 
by : 

Ym(n)  =  cos(ymn) 

Zm(n)  =  sin(ymn)  (5.31) 


where 


w 


m  ~  2 


(m  =  1  ,2, .  . .  ,M) 


(5.32) 


The  above  approximations  assure  that  essential  boundary 
conditions  are  satisfied.  Substitution  of  equations  5.30 
into  5.22  yields: 


KAd  =  AV 


(5.33! 


where 


jT  _  ,,.1.T  /  ,  2 .  T  ,  .m.  T  ,  ,M,T,  .  r 

d  —  1(d)  t  ( d  )  f  • • • f  (d  )  f  • » • f  (d  )  ]  (5.34 


( d  )  =  [u1 ,  v  ^  ^  ^2  9  ^2  9  • • • t  f  « • • r 


ui'  vi] 


(  5.35 


The  stiffness  matrix,  K,  has  the  structure  represented  by 
the  matrix  A,  eq.5.9,  but  in  eq.  5.33  the  submatrix  K1-1  is 
calculated  by 


Klj  =  / (B1) T  6  dn 


(5.36) 


The  vector  of  incremental  pseudo  loads,  for  the  next  time- 
step,  is  calculated  with  the  equation: 


AV  '  =  /B  D  '(A 


T^n  +  1  ,  n+1  .  «n+1 


;HT  tvp  n 


Ax„) dn  +  AFn+1  +  Rn+1 


where  the  residual  out-of-balance  vector  force  is  defined  by: 


„n+ 1  _T  n+1 .  „n+1 

R  =  /B  a  dn  -  F 


(  5.38) 
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6 .  Results  and  Discussions 

6.1-  Reliability  of  the  mathematical  model 

Preliminary  experiments  have  shown  that  the  elasto- 
viscoplastic  response  of  the  FM-73  polymer  layer  can  be 
divided  into  three  main  ranges:  elastic  behavior  at  stress 
levels  below  5  MPa,  primary  creep  -  between  5  and  37  MPa, 
and  secondary  creep  -  above  37  MPa.  The  values  of  the  visco¬ 
plastic  parameters  -  fluidity,  y,  and  power  constant,  N  -  were 
determined  through  curve-fitting  to  the  experimental  results 
obtained  for  the  two  viscoplastic  ranges  and,  as  shown  in 
Figure  26,  good  correlation  was  achieved. 

The  reliability  of  the  numerical  solutions  was  tested 
by  applying  the  present  approach  to  a  wide  class  of  problems 
which  also  have  other  independent  numerical  or  analytical 
solutions.  Here,  the  authors'  main  aim  was  to  find  the 
optimal  number  of  strips  and  global  trial  functions  needed 
to  assure  good  numerical  approximations  to  the  solution  of 
a  wide  variety  of  problems.  As  a  first  step,  the  uncoupled 
diffusion  equation  was  solved  for  the  steady  state  case  with 
constant  coefficients.  The  solution  was  compared  with  well 
known  analytical  solutions  [19].  In  the  case  referred  to  the 
global  trial  functions  in  the  n  direction  coincide  with  the 
eigenfunctions.  Hence,  the  spatial  convergence  in  the  £ 
direction  can  be  examined  by  itself.  Eq.(4.13)  is  of  a 
singular  perturbation  type,  since  the  highest  derivative 
in  the  n  direction  is  multiplied  by  the  small  parameter,  ip . 
Thus,  as  <p  ■*  0 ,  the  interior  solution  approaches  that  of  the 
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one-dimensional  problem  in  the  £  direction.  For  the  time- 
dependent,  one-dimensional  problem  with  nonlinear  coefficients, 
analytical  solutions  are  available  [26], 

Next,  the  convergence  of  the  present  solution  with  that 
of  the  uncoupled  linear  equilibrium  equations,  with  either 
a  uniform  traction  loading  or  a  uniform  temperature  distri¬ 
bution  was  examined.  For  the  cases  of  uniform  loading 
either  mechanical  or  thermal  as  well  as  for  the  case  of  the 
diffusion  equation,  in  which  the  diffusivity  is  function 
only  of  C ,  the  solution  is  composed  of  M  uncoupled  solutions 
i.e.  the  coefficient  matrices  are  of  the  block  diagonal  type. 
Here,  the  effect  of  the  error  introduced  by  the  numerical 
quadrature  can  be  thoroughly  examined.  The  numerical 
solution  of  the  uncoupled  linear  equilibrium  equations,  for 
the  case  of  uniform  temperature  distribution  in  the  NSD 
model,  was  compared  with  the  solution  based  on  the  classical 
theory  of  laminated  plates. 

Finally,  the  convergence  of  the  fully  coupled  nonlinear 
case  was  examined.  The  spatial  discretizations  are  shown  in 
Fig. 27 .  It  was  found  that,  in  uhis  case  (Fig. 28)  as  well  as 
in  the  previously  described  cases,  a  numerical  solution  based 
on  six  strips  gives  accurate  results. 

The  convergence  of  the  time  discretization  was  tested  by 
comparing  results  obtained  after  the  same  total  periods  of 

-3 

time  but  by  different  time  step  increments  (At  =2.10  to 
10_1).  The  solutions  were  found  to  be  unconditionally  stable 
for  0  >  i ,  and  the  maximum  deviation  between  them  was  less 
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than  0.5  per  cent.  The  character  of  the  processes  for  longer 
periods  of  time  is  much  smoother  and  thus  even  more  stable 
and  consistent. 

A  comparison  between  experimental  and  theoretical  results 

for  the  stress  dependency  of  the  dimensionless  diffusivity  is 

shown  in  Fig. 29 .  Experimental  results  were  derived  from 

absorption  tests  under  various  constant  compressive  and  tensile 

stresses  as  described  in  section  3  .  FSM  results  were  obtains-3 

from  successive  solutions,  compressive  or  tensile  stresses 

having  been  imposed  on  the  specimens  at  constant  moisture  and 

temperature.  From  every  solution  the  RRV  was  computed,  while 

the  value  of  the  constant  X_,  for  the  relation,  «  1  +  xdvr< 

o 

was  fitted  only  for  the  linear  relation  zone  (small  stresses) . 
The  good  correlation  shown  in  Fig. 29  implies  that  the  physical 
basis  of  the  RRV  assumption  is  reasonably  solid. 

6 . 2  -Numerical  vs.  experimental  results 

6.2.1;  thick  adhesive  NSD  model 

The  temperature  distribution  along  the  specimen,  at 
several  points  of  time,  is  shown  in  Fig .30  .  a  uniform 
initial  condition  has  been  assumed  for  the  present  problem. 

In  order  to  prevent  the  physically  unrealistic  discon¬ 
tinuities  at  the  free  edges,  initial  conditions  were  smoothed. 
By  using  the  first  four  harmonics,  the  solution  behavior  can 
be  reconstructed.  The  temperature  distribution  across  the 
NSD  specimen's  thickness  at  several  points  of  time  are  shown 
in  Fig. 31.  Since  the  thermal  diffusivity  of  aluminum  is  3 


orders  of  magnitude  greater  than  that  of  FM-73,  the  tem¬ 
perature  distribution  in  the  former  approaches  steady- 
state  more  rapidly  than  in  the  latter.  The  pronounced 
temperature  differences  created  may  thereby  have  a  signifi¬ 
cant  effect  on  the  overall  stresses  and  deformations, 
especially  when  the  visccplastic  behavior  of  the  polymeric 
layer  is  considered. 

The  numerical  prediction  of  the  central  deflection  as 

a  function  of  time,  for  different  values  of  moisture  diffusion 

coefficients,  is  presented  in  Fig. 32.  The  best  correlation 

-5  2 

with  experimental  results  is  obtained  for  D  =  2.5  •  10  mm  /min 

o 

Using  this  value  of  Dq,  several  numerical  solutions  for  the 
central  deflection  vs.  time,  in  two  different  specimens,  are 
shown  in  Fig. 33,  and  these  results  are  compared  with  the 
experimental  findings  described  in  Section  3.  The  FSM 
solutions,  with  and  without  the  moisture-stress  coupling 
effect,  appear  as  the  dashed  and  the  solid  lines,  respect¬ 
ively.  The  coupling  effect  tends  to  improve  the  correlation 
between  numerical  and  experimental  results.  This  improvement 
is  attributable  to  the  fact  that  the  rate  of  diffusion  is 
reduced,  since  the  exposed  polymeric  strip  is  under  compres¬ 
sion,  so  that  the  rates  of  swelling  and  of  straightening  is 
also  reduced.  It  should  be  pointed  out  that  the  direction 
of  the  coupling  effect  changes  in  the  T-1  series  after  about 
fifteen  days.  This  may  be  ascribed  to  the  change  in  the 
direction  of  the  curvature  (Fig. 34),  which  indicates  that 
tensile  stresses  are  prevailing  in  the  outer  polymeric  strips 
(Fig.  35)  . 


moisture  concentration  profiles  are  shown  in  Fin.  35.  The 
average  moisture  concentration  vs.  time  is  shown  in  Fig. 36. 

For  comparison,  experimental  results  are  also  presented. 

The  improvement  in  correlation  due  to  the  use  of  coupled 
effects  is  clearly  demonstrated. 

6.2.2:  thin  adhesive  NSD  model  behavior 

Next,  a  preliminary  investigation  of  3-layers  NSD 
specimens  is  presented.  The  model  is  composed  of  a  0°  araphite/ 
epoxy  layer  which  glues  to  90°  graphite/epoxy  layer  by  means 
of  a  thin  layer  of  FM  73  epoxy  adhesive  (Fig. 37)*.  The  numer¬ 
ical  solution  of  this  case  is  based  on  the  same  procedure 
described  in  section  5. 

The  main  objective  is  to  examine  the  response  of  the 
specimen  due  to  the  viscoplastic  behavior  of  the  FM  73  adhesive. 
Therefore,  for  simplicity,  the  following  assumptions  have  been 
added: 

1 .  Each  lamina  is  described  by  a  macroscopically  homogeneous 
orthotropic  material.  The  material  is  linearly  elastic  under 

a  generalized  Hooke's  law. 

2.  A  2-D  stress  field  is  still  prevailed  in  the  x-y  plane. 

3.  The  problem  is  uncoupled,  namely,  the  equilibrium  equations 
are  solved  for  given  temperature  and  moisture  distributions. 

4.  The  mechanical  loading  is  represented  by  either  a  uniform 
stretching  along  the  y  direction  or  by  a  pure  moment,  along  z 
direction,  achieved  through  a  couple  of  forces  employed  at  the 
"far"  edge. 


The  mechanical  and  hvgro  and  thermo  elastic  parameters  of  the 
graphite/epoxy  and  FM  73  adhesive  age  given  in  Table  6. 
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The  specimen  is  loaded  according  to  the  following  stages: 

1.  A  uniform  thermal  loading  (AT  =  -100°C,  AC  =  0) 

I 

S 

2.  A  given  AC  distribution  through  the  thickness  (see  Fig. 38). 

3.  A  mechanical  loading  which  tends  to  straighten  the 

[  deformed  specimen. 

Numerical  results  of  time  dependency  of  axial  stress  ( a  ) 
distribution  through  the  adhesive  thickness  are  shown  in  Figs. 
39-44.  The  effects  of  uniform  changes  of  temperature,  moisture, 
axial  load  and  moment  on  the  stress  field  vs.  time  functions 
indicate  the  following: 

in  all  cases  relaxation  process  in  the  adhesive  layer  is 
fast  and  stresses  tend  to  relax  to  their  stable  state  in  a  time 

scale  of  less  than  a  minute; 

I 

moisture  tends  to  reduce  significantly  the  tensile  stresses 
in  the  adhesive  even  in  cases  of  mechanical  tensile  loading; 

in  all  cases,  variation  of  axial  stresses  within  the 
adhesive  is  small  and  linear. 

The  present  approach  can  be  directly  applied  to  the 
analysis  of  interlaminar  failures.  In  the  present  problem,  the 
governing  equations  are  of  the  singular  perturbation  type. 

Hence,  for  small  values  of  the  aspect  ratio,  ip ,  it  becomes 
increasingly  difficult  to  solve  these  equations  numerically . 
Several  methods  that  use  singular  perturbation  theories  to 
improve  the  accuracy  and  efficiency  of  the  numerical  solution 
have  recently  been  developed  ([20]  -  [22]).  Incorporation  of 
the  present  approach  into  those  schemes  can  yield  significant 
improvements  for  cases  of  small  aspect  ratios. 


3: 


7 .  Conclusions 

Based  upon  experimental  results  and  numerical  solutions 
of  the  diffusion  and  equilibrium  equations,  the  following  main 
conclusions  can  be  drawn: 

7 . 1  -  The  rate  of  moisture  diffusion  depends  on  the  Relative 

Reduced  Volume  (RRV)  -  the  rate  increases  with  a  rise  in 
the  RRV  due  to  tensile  stress  and  decreases  due  to 
compressive  stress,  and  hygroscopic  swelling. 

7 . 2  -  The  coupling  effect  between  moisture  diffusion  and 

mechanical  behavior  was  determined  quantitatively.  This 
effect  is  pronounced  whenever  the  polymeric  phases  are 
the  major  constituent,  as  in  composites  and  in  thick 
adhesive  layers. 

7 . 3  -  The  mutual  effect  between  moisture,  temperature,  and 

mechanical  behavior,  was  formulated  mathematically. 

The  simultaneous  solution  of  the  nonlinear  coupled 
diffusion  and  equilibrium  equations  provides  an  insight 
into  the  system  which  may  reveal  an  extreme  situation 
in  transition  phenomena. 

7 . 4  -  The  good  correlation  that  was  obtained  between  experi¬ 

mental  and  numerical  results  implies  that  the  proposed 
model  and  approach  are  capable  of  describing  these 
complex  phenomena,  which  are  characterized  by  a  pronounced 
nonlinearity  and  by  time  dependency. 

7 . 5  -  The  finite  strip  method  was  found  to  be  an  accurate  and 

efficient  tool  for  the  solution  of  such  equations. 

7 . 6  -  The  present  study  was  mainly  confined  to  the  nonlinear 

analysis  of  two  material  systems.  However,  it  can  be 
extended  to  the  prediction  of  the  behavior  up  to  failure 
of  any  multimaterial  anisotropic  laminates. 


8 . 1  -  Analytical  study  and  experimental  investigation  of  the 
the  effects  of  hygrothermal  variables  on  failure  and  strength 
of  the  bonded  system. 


8 . 2  -  Investigation  of  the  effect  of  thermal  and  humidity 
fluctuations  (environmental  cycling)  on  stresses  and  strength 
of  the  bonded  system. 

8 . 3  -  Extension  of  the  2-D  analytical  solution  given  in  the 
present  report  into  3-D  analysis  which  is  mainly  applicable 
at  the  free  edges. 

8 . 4  -  Application  of  the  present  approach  and  results  into 
the  analysis  of  interlaminar  behavior  of  composite  laminates. 

9 .  Publications  based  on  present  work 

9.1  -  G.  Yaniv  and  0.  Ishai,  Coupling  Between  Stresses  and 

Moisture  Diffusion  in  Polymeric  Adhesives  (accepted  to 
be  published  in:  Polymer  Engineering  and  Science  Journal 
1986)  . 

9.2 -P.  Bar-Yoseph,  G.  Yaniv  and  0.  Ishai,  The  Interdependenc 
of  Hygrothermal  Processes  and  Elasto-Viscoplastic 
Behavior  in  Polymor-Dominated  Multi-Material  Systers 
(submitted  to  International  J.  of  Solids  and  Structures) 
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stiffness  matrix  components 
strip  lenath  ( FSM  solution) 
flow  vector 

stiffness  matrix  components 

moisture  concentration 

strength  theory  coefficients 

moisture  diffusion  coefficient 

displacement  vector 

elastic  modulus 

effective,  equivalent  stress 

external  tractions  vector 

shear  modulus 

width,  thickness 

stress  invariants 

thermal  conductivity 
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relative  moisture  content  (weight) 
moisture  weight 
local  trial  functions 
compliance  matrix  components 
residual  diffusion  equations 
deviatoric  stress  tensor  components 
relative  reduced  surface 


temperature 


t 


time 


rJ 


physical  variable 

displacement  vector  (x-direction) 

relative  reduced  volume 


u 


v  -  displacement  vector  (y-direction) 
W  -  dissipative  work 

-  weight  of  dry  material 
w  -  displacement  vector  ( z-direction) 
X  -  basis  function  (x-direction) 

Ym  -  global  trial  function 
Zm  -  global  trial  function 


a 

6  - 
Y  " 

i  r  - 
6  - 
£ 

n  - 

9  - 
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thermoelastic  expansion  coefficient 

hygroelastic  expansion  coefficient 

fluidity  coefficient 

region  bound 

central  displacement 

strain  tensor 

non-dimensional  coordinate  (x-direction) 

numerical  summation  parameter 

strain  hardening  coefficient 

compressive-tensile  strength  ratio 

temperature,  moisture,  viscoplastic  work  coupling 
coefficients 

stress-moisture  coupling  coefficient 
m-term  argument 

non-dimensional  coordinate  (in  x-direction) 


^  ,N  .%/•  -V' 
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p  -  density 
o  -  stress  tensor 
t  -  normalized  time 
-  flow  function 

4  -  strip  length  to  width  ratio 

n  -  problem  range 


Abbreviations 

CFRP  -  carbon  fiber  reinforced  plastics 
FEM  -  finite  element  method 
FSM  -  finite  strip  method 
HT  -  hygrothermal 
HOT  -  higher  order  term 
NSD  -  non  symmetrical  doubler 
RMC  -  relative  moisture  content 
RRS  -  relative  reduced  surface 
RRV  -  relative  reduced  volume 
RT  -  room  temperature 
VE  -  viscoelastic 
VP 


viscoplastic 
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Fig.  4  Relative  moisture  concentration  vs.  time  for  epoxy  specimens 

under  different  compressive  loadings 
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Fig.  8  Mechanical  loading  device  for  moisture 
absorption  in  tensile  loading  test 

(see  Fiq.  2)  . 
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different  tensile  loading 
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10  Relative  moisture  concentration  vs.  time  for 
stress-free  FM-73  specimens  with  different 
thicknesses 


Fig.  11  Device  for  measuring  variations  in  central  flexural 

displacement  of  bonded  specimens  during  H.T. exposures 
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Fig.  14  Central  deflection  vs.  time  of  NSD  specimens  (Aluminum  U  FM-73) 
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Fig.  18  Compressive  load  effect  on  relative  moisture  absorption  (Diffusion 
coefficient  variation),  FM-73  specimens  (high  loading) 
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Fig.  20  Tensile  load  effect  on  relative  moisture  absorption 
(Diffusion  coefficient  variation),  FM-73  specimens 
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hygroelastic  compressive  stress  effect  on  relative 
i-n  comparison  to  external  compressive  loads 


Fig.  24  Configuration  and.  boundary  conditions  of  the  domain 

for  the  diffusion  problem 


Fig.  25  Asymmetrical  doubler  -  finite  strip  discretization 
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Fig.  29  Dependence  of  relative  diffusivity  on  stress  (1st  invariant  J 
comparison  between  expr.  results  and  numerical  evaluations  of 
the  influence  of  RRV  variations 
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Fig.  35  Schematic  description  of  moisture  and  stress  thickness 
distributions  and  deflections  along  the  specimen  at 
different  time  levels  (  T=-39°C) 
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Fig.  37  A3  layers  thin  adhesive  NSD  specimen  geometry 

and  properties 


Fig.  38  Typical  moisture  profiles 
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Moisture  distribution  over  specimen  width  at  free  absorption 
4yyroelastic  compression  zone 


Table  1 ;  Hygrothermomechanical  Conditions  for  Epoxy  Specimens 
(Compression) 


Comments 

Test  Duration 
(no.  days) 

Compressive 

Load 

(MPa) 

Test 

Temperature 

(°C) 

Series 

3  specimens 

30 

0 

40 

E-1 

3  specimens 

30 

1  .8 

40 

E-2 

2  specimens 

30 

3.5 

40 

E-3 

3  specimens 

30 

5.2 

40 

E-4 

Failure 

14 

7.6 

40 

E-5 

Table  2:  Hygrothermomechanical  Conditions  for  FM-73  Specimens 
Compression  (Low  Loading) . 


Comments 

Test  Duration 
(no.  days) 

Compressive 

Load 

(MPa) 

Test 

Temperature 

(°C) 

Series 

4  specimens 

1  4 

0 

40 

FM-  1 

4  specimens 

14 

2 

40 

FM-2 

4  specimens 

14 

4  • 

40 

FM-3 

Failure 

4 

6 

40 

FM-4 
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Table  3 ;  Hygrothermomechanical  Conditions  for  FM-73  Specimens 
Compression  (High  Loading). 


Comments 

Test  Duration 
(no.  days) 

Compressive 

Load 

(MPa) 

Test 

Temperature 

(°C) 

Series 

2  specimens 

60 

0 

40 

H-1 

2  specimens 

60 

8 

40 

H-2 

2  specimens 

60 

155 

40 

H-3 

2  specimens 

8 

20 

40 

H-4 

Failure 

2 

25 

40 

H-5 

Table  4;  Hygrothermomechanical  Conditions  for  FM-73  Specimens 
(Tension) 


Comments 

Test  Duration 
(no.  days) 

Tensile 

Load 

(MPa) 

Test 

Temperature 

(°C) 

Series 

4  specimens 

21 

0 

40 

TM-1 

4  specimens 

21 

4 

40 

TM-2 

4  specimens 

21 

8 

40 

TM-3 

4  specimens 

21 

12  . 

40 

TM-4 

4  specimens 

21 

1  6 

40 

TM-5 

Failure 

6 

20 

40 

T-6 

Order 

Adherend 

1 

Adherend  2 

T 

SF 

[°C] 

Material 

B 

Material 

mm 

T-1 

At-2024  T3 

1 

FM-73 

3.5 

B 

C-1 

A£-2024  T3 

1 

FM-73 

3.5 

C-2 

AA-2024  T3 

1 

FM-73 

5.2 

B 

Table  5 


Material  composition  and  dimensions 
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Appendix  A  -  Method  for  Determination  of  Relative  Moisture  (RMC) 
Absorbed  by  Specimens  Under  Tension 

The  RMC  was  calculated  as  follows: 

The  weight  of  a  dry  specimen  is  given  by: 

Wd  =  EWd  +  IWd  (A- 1 ) 

where  £ W d  and  are  the  dry  weiahts  of  the  end  and  inner 

zones,  respectively. 

The  total  weight  of  a  wet  specimen  is 

W  =  EW  +  ];W  (A- 2) 

and  the  absorbed  moisture  is: 


m  =  £m  +  Im  .  (A-3) 

The  RMC,  calculated  for  the  uniformly-stressed  inner  zone,  is 
given  as: 


M  = 


IWd 


or,  for  the  k-th  specimen  (k  =  1,2, 3, 4) 


M  = 


k  k 
m  -Em 

W*5 

IWd 


(A-4) 


(A- 5) 


Since  the  end  zones  are  not  included,  the  RMC  of  the  k-th  and 
the  k+1  specimens  (in  the  same  test)  should  be  the  same: 


Mk  =  Mk+1 


( A-6 ) 


Substituting  eq.  (A-5)  in  (A-6) ,  and  solving  for  £m  we  obtain: 

(A- 7) 


where 


.rr 


k  +  1 


/a"  \ '  •.* 


*  c.* 


*  *  *  » ^  «  *  i  *  »  '  «  »  >  ’*■ 
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Substitution  of  m  from  eq.  (A-7)  in  eq.  (A-5)  yield: 


Mk  = 


IW**1  .  -  VI 

_ Q _ a__ 

<wd+1  -  wdJ 


( A-8) 


or: 


_AW 

AW 


(A-9) 


where  AW  is  the  weight  difference  between  two  wet  specimens 
in  the  same  test,  and  AW^  the  dry  weight  difference  between 
them.  It  should  be  noted  that  the  result  for  each  test  is  an 
average  of  six  values,  although  there  are  four  specimens. 
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Appendix  B  -  Evaluation  of  Equivalent  Hvaroelastic  Stress 


in  Unstressed  Thick  Specimens 

Let  the  distribution  of  moisture  concentration  at  a 
certain  time,  C(y,t),  be  described  by  Fig. 45a.  The  dashed 
line  designates  the  average  concentration,  C.  The  point  SFC 
(stress-free  concentration)  divides  the  thickness  of  the 
specimen  into  a  hygroelastically  compressed  region  - 
y  ^  y  S  h/2  with  moisture  concentration  above  C,  and  a  tensile 
region  0  S  y  i  y  with  concentration  below  C. 

The  total  axial  strain  of  a  layer  dy  at  distance  y  from 
the  centerline  of  the  specimen  is  given  by: 

ex(Y)  =  ee  +  eH  (B~1) 

or 

ex(y)  =  |Ox  (y)  +  SC  (y)  (B-2) 

v/here  is  the  axial  normal  stress,  B  the  coefficient  of 
hygroelastic  expansion  (CHE) ,  and  C  is  moisture  concentration. 

The  compressive  stress  required  to  reduce  the  strain, 
within  the  range  y  £  y  i  h/2  to  its  average  level,  e  is: 

cx(y)  =  E [ SC ( y )  -  e]  ,  (B-3) 

or,  with  the  average  strain  expressed  in  terms  of  the  average 

concentration  C: 

ox(y)  =  SE[C(y)  -  C]  .  (B-4) 


Let  us  concentrate  on  the  compressed  zone,  Fig. 4 5b,  and 


define  a  new  coordinate  system  as  follows: 

„  h 

*  =  2  "  y  (B-5) 

C*  (Y)  =  C  (y )  -  c 


According  to  equations  B-5,  the  equivalent  compressive  stress 
can  be  defined  as: 


=  -  /  0EC* (Y) dY 
Y  0 


(B-6) 


Assuming  0  and  E  to  be  independent  of  Y,  we  may  approximate: 


a  =  48EC* 
x  j  rn 


(B-7 ) 


Transformation  to  the  original  coordinate  system  yields  the 
expression  for  the  equivalent  stress: 
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Appendix  C  -  The  Procedure  of  the  Computer  Program 

The  program  relates  to  a  structure  with  rectangular 
bounds  composed  of  general  layers  of  different  materials, 
which  is  exposed  to  a  given  hygrothermomechanical  external 
conditions.  The  program  is  able  to  compute  the  tempterature , 
moisture  and  mechanical  fields  which  are  developed  within  the 
layered  structure  and  the  interaction  -  coupling  between  all 
these  process  variables.  The  computer  program  is  composed  of 
two  main  parts  -  the  first  solves  the  diffusion  equations 
whereas  the  second  computes  the  mechanical  response  cf  the 
structure  based  on  assumption  of  visco-elasto-plastic 
behavioral  model.  The  first  part  is  solved  iteratively 
up  to  the  level  of  desired  convergence,  (Modified  Newton- 
Raphson ,  Newton-Raphson-Picard) .  In  the  second  part,  the 
residual  load  (or  the  out-of-balance  load)  is  added  to  the 
load  vector  of  the  next  time  step. 

The  numerical  solution  is  conducted  by  the  finite  strip 
method  (FSM) .  The  division  of  the  structure  into  longitudi¬ 
nal  strips  is  done  automatically,  whereas  there  is  a  possi¬ 
bility  of  obtaining  three  different  division  regions  through 
the  width  of  the  structure.  The  proportions  between  the 
different  division  regions  and  the  number  of  elements  in  each 
region  are  free  to  be  determined  according  to  the  given 
requirements.  In-case  of  bi-material  structure,  one  region 
at  each  side  of  the  interface  is  selected.  To  improve  the 
precision  of  the  solution,  the  strips  were  divided  into  two 
sub-regions  at  which  an  additional  division  to  sub-strips  of  a 
different  number  and  density  is  possible. 


TF&mrsm . . . 
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In  the  present  work,  the  integrals  were  calculated 
numerically  according  to  the  Gauss-Legendre  scheme.  This  was 
done  for  three  integration  points  in  each  direction  through 
the  strip  width  and  along  the  sub-strip,  to  obtain  optimal 
adjustment  for  computation  for  square  interpolation  functions 

The  computer  program  starts  with  the  determination  of 
dynamic  dimensions  for  most  of  the  main  systems  to  save 
memory  resources  (use  of  memory  volume  is  required  without 
useless  surplus) .  An  additional  means  of  saving  memory 
resources  is  by  determination  of  joint  working  spaces. 
Detailed  illustration  of  the  computer  program  scheme  and 
steps  is  given  in  the  following  flow  charts  (Fig. 46). 
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